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ABSTRACT 


How diffuse dwarf galaxies that are deficient in dark matter—such as NGC1052-DF2 and NGC1052- 
DF4—formed is a mystery. Along with their luminous member globular clusters (GCs), the so-called dark 
matter deficient galaxies (DMDGs) have challenged observers and theorists alike. Here we report a suite of 
galaxy collision simulations using the adaptive mesh refinement code ENZO with 1.25 pc resolution, which 
demonstrates that high-velocity galaxy collisions induce the formation of DMDGs and their star clusters (SCs) 
simultaneously. With a numerical resolution that is significantly better than our previous study, we resolve the 
dynamical structure of the produced DMDGs and the detailed formation history of their SCs, which are possible 
progenitors of the DMDG's member GCs. In particular, we show that a galaxy collision with a high relative 
velocity of ~ 300 km s~!, invoking severe shock compression, spawns multiple massive SCs (M, > 106 Mo) 
in « 150 Myr after the collision. At the end of the ~ 800 Myr evolution in our fiducial run, the resulting 
DMDG of M, ~ 3.5 x 108 Mo hosts 10 luminous (My X —8.5 mag), gravitationally bound SCs with a line-of- 
sight velocity dispersion 11.2 km sl. Our study suggests that DMDGs and their luminous member SCs could 
form simultaneously in high-velocity galaxy collisions while being in line with the key observed properties of 
NGC1052-DF2 and NGC1052-DF4. 
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1. INTRODUCTION 


At the center of the intense debate on the formation mech- 
anism that can address abnormal lack of dark matter and un- 
usually bright populations of globular clusters (GCs) are two 
recently discovered ultra-diffuse dwarf galaxies, NGC1052- 
DF2 and NGC1052-DF4 (hereafter DF2 and DF4, respec- 
tively). The spectroscopic study of their member GCs has 
revealed surprisingly small velocity dispersions compared to 
what is normally expected by their observed stellar masses. 
This implies that DF2 and DF4 have very small dark-matter- 
to-stellar-mass ratios (dark matter deficient galaxy (DMDG); 
van Dokkum et al. 2018b, 2019).! The exceptionally large 
population of their luminous member GCs is also a mys- 
tery. The GC luminosity function (GCLF) of DF2 peaks at 
My ~ —9, approximately 1.5 mag shifted to the brighter side 
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from the Milky Way's GCLF (van Dokkum et al. 2018c; Ma 
et al. 2020a; Shen et al. 20212), indicating a large number of 
massive GCs. Extensive spectroscopic studies of DF2 have 
found that their stellar population is ~9 Gyr old on average 
and that the metallicity of the main stellar body differs from 
that of its member GCs (Fensch et al. 2019; Ruiz-Lara et al. 
2019). As for the kinematics, recent studies hint that DF2 
is likely a prolate rotator (Emsellem et al. 2019; Lewis et al. 
2020). 

These observations could provide a testbed for various for- 
mation scenarios proposed for DF2 and DF4. In one of such 
scenario suggested by Silk (2019), a “mini-Bullet cluster- 
like" event on a galaxy scale decouples dark matter from 
dissipative baryons.” The high-velocity galaxy collision in- 
duces a strong shear and shock compression in the interstel- 
lar medium (ISM), instigating the formation of DMDGs and 
their bright member GCs (Silk 2019; Trujillo-Gomez et al. 
2020). Numerical studies have demonstrated the formation 


1 This anomalous characteristic is based on the distance measurement, the 
accuracy of which is still under debate (see e.g., Trujillo et al. 2019 and 
Monelli & Trujillo 2019 for details, but also van Dokkum et al. 2018a, 
Danieli et al. 2020 and Shen et al. 2021b). 


? Other scenarios include the tidal dwarf galaxies (Fensch et al. 2019) and the 
dwarf galaxies with tidally stripped dark matter halos (Ogiya 2018; Jackson 
et al. 2021). See Shin et al. (2020) for more information. 
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Table 1. A Suite of Idealized Galaxy Gollision simulations listed with their initial configurations and the stellar mass of the most massive DMDG formed 


Run Name Stellar Feedback Energy ^ Initial Gas Metallicity Disk Inclination Angle Disk Spin Relative Angle M, pMDG tend 
(10°! erg/150 Mc, of Stars) (Zo) (81,01), (62.02) (108 Mo) (Myr) 
a) (2) (3) (4) (5) (6) (7) 
Fiducial/ 1 1 (0°, 0°), (0°, 0°) 0° 3.47 (3.72) 800 (250) 
FB3 3.3 1 (0°, 0°), (0°, 0°) 0? 2.76 250 
FB7 6.7 1 (0°, 0°), (0°, 0°) 0° 1.46 250 
FB10 10 1 (0°, 0°), (0°, 0°) 0° - 50 
Low-metal* 1 107 (0°, 0°), (0°, 0°) 0° 0.90 250 
Retrograde 1 1 (0°, 0°), (180°, 0°) 180° 1.76 250 
Tilt-1 1 1 (—22°5, 0°), (22°5, 0°) 45° 2.33 250 
Tilt-1-FB10 10 1 (—22°5, 0°), (22°5, 0°) 45° - 50 
Tilt-2 1 1 (0°, —22°5), (0°, 22°5) 45° - 50 
Tilt-2-FB10 10 1 (0°, —22°5), (0°, 22°5) 45° - 50 


NOTE—The simulation configurations and the stellar mass of the most massive DMDG formed: (1) run name, (2) stellar feedback energy, (3) initial metallicity of the gas disks, (4) 


angle by which each of the two colliding galactic disks is rotated about the +x-axis (01, 03) and about the +y-axis (91, $2), starting from its original position in the x-y plane (see the 


top-left panel of Figure 1), (5) relative angle between the angular momentum vectors of the two galactic disks, (6) stellar mass within 12 kpc from the center of the most massive 


DMDG (“-” = no DMDG formed), (7) time since the pericentric approach of the two progenitor disks when we end the simulation. 


+ This is the fiducial run presented in all of the figures in this article. See Sections 2 and 3 for more information. 


$ This is the run presented in the right panels of Figure 5 in Section 3.2. 


of young massive clusters (YMCs) and/or GC candidates 
triggered by galaxy mergers in idealized simulations (Li et al. 
2004; Bournaud et al. 2008; Kruijssen et al. 2011, 2012; Hop- 
kins et al. 2013; Renaud et al. 2015; Maji et al. 2017; Moreno 
et al. 2019; Lahén et al. 2019, 2020) and in cosmological sim- 
ulations (Li et al. 2017; Kim et al. 2018; Ma et al. 2020b). 
To validate the scenario by Silk (2019), we previously simu- 
lated idealized galaxy collisions using two different gravito- 
hydrodynamics codes and showed that ~108 Mo DMDGs are 
produced when two gas-rich, dwarf galaxies collide with a 
relative velocity of ~ 300 km s7! (Shin et al. 2020). How- 
ever, the chosen resolution of 80 pc limited our ability to 
examine the inner structure of the resulting DMDGs, or to 
resolve the formation of their member SCs. 

In this Letter, we examine the galaxy collision-induced 
DMDG formation scenario using adaptive mesh refinement 
(AMR) simulations with significantly improved 1.25 pc res- 
olution. We not only demonstrate the DMDGs' formation but 
also explore the detailed formation histories and the dynam- 
ical properties of the DMDG and their member SCs which 
could eventually evolve into bright GCs. 


2. SIMULATIONS 


The initial condition used for this study is similar to what 
was adopted in Shin et al. (2020) but with greatly improved 
particle resolution: a dark matter halo of mass 3.88 x 10? Mo 
with a particle mass 1.55 x 10? Mọ each, a stellar disk of 
4.21 x 105 Mo with a particle mass 1.05 x 10? M. each, and 
a gas disk of 1.68 x 10? M. The initial gas fraction is 
thus fas = Mgas/ Mgaiaxy = 0.28, but because ~ 30% of the 
gas is consumed before the collision, fgas becomes ~ 0.2 
when two galaxies collide at t = 0. Our Fiducial run is 
a co-planar prograde-prograde collision of the two identical 


galaxies with a relative velocity of ~ 300 km s at a 40 kpc 
distance, and with a pericentric distance of 1 kpc (see Ta- 
ble 1)? By varying the stellar feedback energy, initial gas 
metallicity, disk inclinations and spins, we form a suite of 
numerical experiments. The simulation configurations and 
the stellar masses of the resulting DMDGs are listed in Table 
jm 

The simulations are performed with the publicly available 
code ENZO (Bryan et al. 2014; Brummel-Smith et al. 2019) 
with the ZEUS hydrodynamics solver (Stone & Norman 
1992a,b). Radiative gas cooling and heating rates along with 
the metagalactic ultraviolet (UV) background are computed 
by the GRACKLE library that interpolates the lookup table 
from CLOUDY (Smith et al. 2017). To describe feedback- 
regulated star formation, we employ the star-forming molec- 
ular cloud model (SFMC; for details, see Kim et al. 2013, 
2019), a change from a simpler prescription in Shin et al. 
(2020). A particle representing a SFMC forms when (1) the 
density of a gas cell exceeds 4 x 10^ cm ^, (2) the gas flow 
is converging, (3) the cooling time of the cell is shorter than 
its dynamical time (f4yn), and (4) the created SFMC particle 
is heavier than Msrmc = 10? Mo. Depicting inefficient star 
formation in molecular clouds (Krumholz & Tan 2007), the 
SFMC particle returns 8096 of its original mass in 12 lbi" 
together with thermal energy peaking at I tayn (four different 


? Readers should note that this initial condition models a rare type of hy- 
perbolic galaxy collision events, but a type that does exist in a simulated 
universe such as IllustrisTNG (for details, see Shin et al. 2020). 

^ The local star formation efficiency (SFE) per freefall time (£r ~ 1%) is 
based on Kim et al. (2011, 2019). Other studies have also found that vary- 
ing the SFE parameter by not more than an order of magnitude does not 
significantly affect the SC properties (Li et al. 2020; Ma et al. 2020b). 
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Figure 1. ¢ = —100, 50, 150, 800 Myr snapshots of the fiducial, idealized collision of two identical gas-rich dwarf galaxies of M200 = 


5.95 x 10? Mo each, and the resulting DMDGs (the Fiducial run in Table 1). Surface densities of gas (top row), stars (middle row) and dark 
matter (bottom row). t = 0 is set to the moment when the two progenitor galaxies are at pericentric approach. The co-planar galaxy collision 
(i.e., both progenitor disks lying in the x-y plane) with a relative velocity of 300 km s produces multiple DMDGs with stellar mass ~108 Mo 
each. The dotted circles of 10 kpc radius are centered on the most massive DMDG produced. See Section 3.1 for more information. 


energy choices for Fiducial to FB10 runs between 10?! 3? 
erg per 150Mọo of stars formed) and metals (2% of the re- 
turned gas; see also Cen & Ostriker 1992; Kim et al. 201 1).5 


A pair of colliding galaxies is initialized in a (1.311 Mpc)? 
box of 64? root grid. We adaptively refine the galaxies up to 
level /max = 14 achieving 1.25 pc resolution.* After the two 
progenitor galaxies pass each other (¢ = 100 Myr), we re- 


5 While it is possible that the momentum imparted by supernovae could be 
underestimated in our thermal feedback model, the subgrid physics imple- 
mentation adopted here have consistently shown that the the star forma- 
tion is efficiently self-regulated and the DMDGs are formed, across multi- 
ple resolution (resolution study at 80, 5, 1.25, and 0.625 pc) and feedback 
strengths (Fiducial, FB3, and FB7 run in Table 1). 


6 We create a relaxed galactic disk by gradually increasing the resolution for 
250 Myr, thereby suppressing artificial starbursts: maximum spatial reso- 
lution = 80 pc for the first 100 Myr (/max = 8), 20 pe for the next 100 Myr 
(Imax = 10), and 5 pc for the last 50 Myr (Imax = 12). The relaxation step 
ends at t = —150 Myr, the beginning of our 1.25 pc resolution simulation 
discussed in this article (t = 0 is when the two galaxies are at pericentric 
approach; see, e.g., Figure 4). During the relaxation, SFMC formation is 
either not allowed (i.e., dense gas clumps are supported only by the Jeans 
pressure floor) or allowed with a coarser-resolution recipe. 


fine the cells only in a pre-defined region that only includes 
the produced DMDGs. This RefineRegion is initially set 
to a (16 kpc)? box containing the most massive DMDG 
throughout the run, and is surrounded by static nested boxes 
of cascading resolution from level 0 to 8. At level /, a 
cell splits into eight child cells whenever the cell contains 
more gas than MZ, gas —2 05540 1) x Mid as 
ticles than Ml. part = ge x p14 , Where M 


ref, pari ref, gas T 
2000 Mo = 2 MER and Mt pan = 4000 Mo = 4 Msrwc- 


or more par- 


3. RESULTS 


3.1. Collision-induced DMDGs: Kinematics and 
Morphology 


Figure 1 shows the collision sequence of two gas-rich 
dwarf galaxies approaching with a relative velocity of ~ 
300 km s7! (the Fiducial run in the first row of Table 1; 
similar to what was extensively tested in Shin et al. 2020). 
We define t = 0 as the time when the two galaxies are at 
pericentric approach. As shown in the top row, we find that 
extremely dense gas clumps form during t = 0 — 50 Myr as 
the colliding columns of gas undergo severe shock compres- 
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Figure 2. Left panel: the radial profile of enclosed masses in a logarithmic scale from the center of the most massive DMDG at t = 800 Myr 
in the Fiducial run. The DMDG consists of mostly stars, lacking in dark matter and gas. The vertical magenta line denotes the half-mass 
radius Rh of the DMDG. Middle panel: the dynamical mass within r estimated from the individual SC's circular velocity, M< (r) = rv2,,.(r)/G, 
now in a linear scale plot (circles). Among 42 gravitationally bound SCs the DMDG hosts, the five most massive SCs are shown with colored 
markers. For comparison, the true curve of the enclosed total mass is also included (dotted line with squares). Right panel: the line-of-sight 
stellar velocities in a (5 kpo)? box centered on the DMDG at t = 800 Myr, projected to a plane whose vertical coordinate (z^) is aligned with 
the total angular momentum of the DMDG (Lpmpo). Colored circles denote the location of SCs and their line-of-sight velocities with radii of 


2 (Msc/ Mo)!’ kpc. Two isosurface density contours are also shown: log[E, / (Mz; pc ?)] = 0 and 1. See Section 3.1 for more information. 


sion near the first contact point (the center of each panel) 
and along the tidally stripped gas streams. The highly com- 
pressed gas clumps of X4, ~ 10?-^ Mc pc? (see also Fig- 
ure 3)—a typical value for extremely dense giant molecular 
clouds—invokes the formation of DMDGs and their mem- 
ber SCs (Kim et al. 2018; Lahén et al. 2020; Ma et al. 
2020b). Strong bursts of star formation continue inside the 
gas clumps during t = 50 — 150 Myr (middle row). The newly 
formed stars remain near the first contact point and begin to 
form one or more gravitationally bound structures without 
the help of pre-existing dark matter potentials. In contrast, 
the two dark matter halos simply pass through each other dur- 
ing this high-velocity collision. As a result, as seen in the bot- 
tom row of Figure 1, hardly any dark matter particle remains 
in the resulting DMDGs. Among the three DMDGs identi- 
fied by the Hop halo finder (Efstathiou et al. 1985) in the 
Fiducial run, the most massive one has M, = 3.5 x 108 Mo 
and almost no dark matter (see also Figure 2). In Figure 1, 
the most massive DMDG found is marked with a dotted cir- 
cle in each panel. The DMDG drifts away slightly from the 
first contact point by t = 800 Myr (the fourth panels in each 
TOW). 

In Figure 2, we present the dynamical properties of the 
most massive DMDG produced in the Fiducial run at 
t — 800 Myr. The left panel displays the radial profiles 
of the enclosed masses in different components: stars, gas, 
and dark matter. The DMDG primarily consists of stars 
and is nearly devoid of not only dark matter but also gas, 
which is in line with the reported gas-poor nature of DF2 
(Mas, DF2 < 3.15 x 109 Mo; Chowdhury 2019). The half- 


mass radius Åh of the DMDG's stellar body is 0.026 kpc, 
which is > 10 times smaller than that of a typical ultra-diffuse 
galaxy (~ 1.5 kpc) or that of DF2 (~ 2.2 kpc; van Dokkum 
et al. 2015, 2018b). However, it is important to note that our 
simulated DMDGs evolved only for « 800 Myr in an ideal- 
ized, isolated environment. Some morphological properties 
of a compact stellar object, such as Mas and Rp could be sen- 
sitive to the energetic supernovae feedback (Trujillo-Gomez 
et al. 2021), or the tidal interaction with massive neighbors." 
Rh also changes frequently at this early stage of evolution as 
a number of new SCs in the outskirts leave the DMDG, or 
because of the dynamical relaxation or heating. 

We now study the kinematics of SCs formed simultane- 
ously with the DMDG that are identified by the HOP algo- 
rithm reconfigured to locate SCs. We fit each SC's stellar 
density profile to a power law, p(r) = po (14-r2/a2) (Y*0? 
(Elson et al. 1987), and obtain its scale radius a. R — 10a is 
adopted as the extent of each SC when we determine its phys- 
ical quantities such as mass, half-mass radius, orbital veloc- 
ity, average metallicity, and creation time. In the Fiducial 
run, a total of 50 SCs are found within 12 kpc from the 
center of the DMDG at t = 800 Myr; among them, 42 are 
gravitationally bound to the DMDG, and eight are unbound. 
The middle panel of Figure 2 presents the dynamical mass 
within radius r estimated from each SC's circular velocity, 


1 For example, the DMDGs formed in an isolated simulation may harbor 
more gas than DF2 does. However, Shin et al. (2020) demonstrated that 
if the simulated DMDG orbits around another massive galaxy, the tidal 
interaction could deprive the DMDG of most of its gas in a few Gyr. 
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Figure 3. Similar to Figure 1, but zoomed in on the most massive SC (*SC1”) created around the most massive DMDG in the Fiducial run. 
t = 70, 110, 160, and 210 Myr snapshots of the gas surface density and stellar distribution in a (1 kpc)? box (top and middle rows), and in a (6 
kpc)? box (bottom row). The red dots in the middle and bottom panels denote the stars in SC1 (the same color codes and SCs as in Figure 2). 
The dotted circle of I kpc radius marks the location of the DMDG in each of the bottom panels. See Section 3.2 for more information. 


M. (r) =rv2.,.(r)/G. Shown for comparison is the true curve 
of the DMDG’s enclosed mass at r. It is notable that the dy- 
namical masses estimated by the gravitationally bound SCs 
lie close to the true enclosed mass. This implies that the grav- 
itationally bound “member” SCs of the DMDG are virialized, 
a key assumption when the dynamical masses of stellar sys- 
tems are observationally evaluated. 

In the right panel of Figure 2, we inspect the dynamical 
structure of the DMDG in the Fiducial run. It displays the 
line-of-sight velocities of the stars and SCs in the DMDG, 
viewed from an angle that is perpendicular to the angular 
momentum vector of the DMDG (Lpwpo) at t = 800 Myr. 
We first find that the line-of-sight velocity dispersion of the 
42 gravitationally bound SCs is o = 12.6 kms^!. The 10 
most luminous bound SCs (My < —8.5; see Figure 6 and 
Section 3.2) have o = 11.2 km s™!. The small spread in ve- 
locity—smaller than what is normally expected by its stellar 
mass—is an indication that the resulting galaxy has a small 
dark-matter-to-stellar-mass ratio. Our finding is in line with 
the observed value of 8.4 km s^! for the luminous SCs in 
DF2 (van Dokkum et al. 2018b) and 5.8 km s lin DF4 (van 
Dokkum et al. 2019), considering that the observation is con- 


ducted in an oblique plane. Second, from the “butterfly- 
shaped" line-of-sight velocities and the isosurface density 
contours, we find that the DMDG is an oblate rotator with 
a major-to-minor axis ratio ~ 2. The rotation velocity of the 
stars is ~ 15 kms ^l at 4 Rh < r < 10 Rh with a noticeable ve- 
locity gradient along the major axis of the system at r < Rh 
(Veire,max ^v 35 km s latr-—25 pc). For comparison, the 
stellar body of DF2 is reported to be a prolate rotator with 
a rotational velocity of 10.8 — 12.4 km s^! (Emsellem et al. 
2019 and Lewis et al. 2020, respectively). 

However, by comparing the specific angular momentum 
of the DMDG (/pwpo) and that of the system of 42 bound 
“member” SCs around the DMDG's center (sc), we dis- 
cover that the SCs are not in a coherent rotation around the 
DMDG, nor preferentially in the disk plane of the DMDG. 
(i.e., |Isc|/ Ipupa| = 0.042). The rotation axis of the system 
of SCs about the DMDG—f it exists at all—is not aligned 
with that of the whole DMDG (i.e., Isc and Ipupg makes 
a 69? angle). The discrepancy in the dynamical states of the 
DMDG and the system of its member SCs could be attributed 
to their different formation processes (see Section 3.2). 
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Lastly, using the suite of simulations that we performed 
(Table 1), we remark on how the formation of the collision- 
induced DMDGs depends on the simulation parameters. 
(1) First, DMDGs form if a sufficiently large amount of 
gas participates in the supersonic collision and forms self- 
gravitating gas clumps. For example, in the Tilt-1 run, 
although the collision is no longer co-planar, a portion of 
the disk gas along the x-axis (see Figure 1) still clashes 
into a large column of gas in the other disk, producing two 
DMDGs. In contrast, in the Tilt-2 run, the disk gas finds 
only a small column of gas in the other disk during its first 
near head-on collision, so no DMDG forms (for a quanti- 
tative discussion, see Section 3.3 of Shin et al. 2020). (2) 
The strength of stellar feedback affects the internal and ki- 
netic energy of the ISM and dictates the DMDG formation. 
In the FB10 run where the feedback is 10 times more ener- 
getic than in the Fiducial run, clumps are easily dispersed 
via stellar feedback, with no DMDG formed. (3) The metal- 
licity of the disk gas influences the DMDG formation. In 
the Low-metal run where the initial gas metallicity is set to 
107? Zs (100 times lower than in the Fiducial run), the rel- 
atively less metal-enriched gas slows down the gas cooling 
and suppresses star formation. As a result, the mass of the 
resulting DMDG is only one-third of what is found in the 
Fiducial run, hosting only nine SCs (as opposed to 50 in 
the Fiducial run; see also Figure 5). (4) The velocity shear 
in the gas at the first contact surface affects the dynamical 
properties of the resulting DMDG. In the Fiducial run of 
a prograde-prograde collision, the ISM experiences a strong 
shear as the colliding columns of gas move in the opposite di- 
rections. Therefore, a DMDG is born in a clump that is in co- 
herent rotation, resulting in a rotation-supported structure. In 
contrast, in the Retrograde run with a prograde—retrograde 
collision, the velocity vectors in the collapsing gas clumps 
are highly disordered, resulting in a DMDG with less rota- 
tion (i.e., |lRetrograde|/|lFiducia1] i 0.1). 


3.2. SCs in the Collision-induced DMDGs: Formation 
History and Key Observables 


We turn our attention to the formation history and the prop- 
erties of SCs in the collision-induced DMDGs. Figure 3 
again shows the collision sequence of the Fiducial run, but 
zooming in on the site where the most massive SC of the 
most massive DMDG was formed (“SC1”; red dots). After a 
large fraction of the disk gas loses its initial momentum dur- 
ing the first near head-on collision and stalls at the first con- 
tact surface, pockets of dense, shock-compressed gas emerge 
(t = 50 Myr panels in Figure 1). A few of these very dense 
clumps (Zgas Z 10?-^ M; pc?) rapidly cools and collapses 
to spawn clusters of stars before or while the DMDG's main 
stellar body forms (M, sc Z 109 Mc each; t = 70 Myr pan- 
els in Figure 3; Kim et al. 2018; Madau et al. 2020). As the 
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Figure 4. Top panel: the star formation history in the most massive 
DMDG and its five most massive member SCs in the Fiducial 
run (the same SCs and colored markers as in Figure 2). Bottom 
panel: zoomed in on the time that the five SCs are created. The 
five SCs and most of the stars in the DMDG's stellar body form 
during t = 0 — 150 Myr, immediately after the progenitor galaxies’ 
pericentric approach. See Section 3.2 for more information. 


clumps continue to turn gas into stars in the next ~ 50 Myr, 
they are fed by the accreting cold gas while the frequent su- 
pernova explosions from young SFMCs eject the hot gas out 
of the new SCs (t = 110 Myr panels in Figure 3; see also the 
bottom panel in Figure 4). The SCs consume or lose most 
of their gas by t = 160 Myr, and continue to orbit around the 
DMDG until the end of our simulation at t = 800 Myr. 

Figure 4 presents the star formation history in the most 
massive DMDG of the Fiducial run and its five most mas- 
sive member SCs. Most of the stars in the DMDG and the 
five SCs are born in < 150 Myr after the progenitor galaxies’ 
pericentric approach (t = 0-150 Myr), which is consistent 
with our earlier studies (Shin et al. 2020). After the initial 
burst of star formation, the DMDG maintains to form stars at 
a modest pace (~ 107! — 107? Mo yt"), while none of the 
massive SCs form additional stars at a meaningful rate after 
t = 150 Myr as they are nearly devoid of gas. The widths of 
star formation peaks for these SCs are ~50 Myr, which are 
significantly narrower than that of the DMDG. 

Then in Figure 5, we make scatter plots of all the SCs 
formed with the DMDG on the plane of the stellar mass ver- 
sus the average metallicity, and on the plane of the average 
creation time versus the average metallicity. The left pan- 
els are for the Fiducial run and the right panels are for 
the Low-metal run (note that the Low-metal run with the 
initial metallicity of 107? Zo produces nine SCs compared 
to 50 in the Fiducial run; see Table | and Section 3.1 for 
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Figure 5. Top row: the stellar masses versus the average metal- 
licities of the SCs formed in the DMDG. Bottom row: the average 
creation times versus the average metallicities of the SCs. The left 
panels are for the Fiducial run (42 bound 4- eight unbound SCs) 
and the right panels are for the Low-metal run (two bound + seven 
unbound SCs). Selected massive SCs are shown with the colored 
markers (the same SCs and markers as in Figures 2 and 4 for the 
Fiducial run) The average metallicities of SCs are correlated 
with their masses. We also denote the average creation time and 
metallicity of the DMDG as a whole with a black circle in the bot- 
tom panels. In the Low-metal run, the stars in the DMDG is more 
metal-enriched than the two bound SCs. See Section 3.2 for more 
information. 


more information). In the bottom panels, we also denote the 
average creation time and metallicity of the stellar body of 
the DMDG with a black circle. Two observations could be 
made. First, we find that the average metallicity of a SC cor- 
relates positively with its mass (top row) rather than with its 
creation time (bottom row). We argue that because a massive 
SC creates more stars that consequently injects more met- 
als into the ISM, and retains a higher fraction of supernova 
ejecta in its deeper gravitational potential well, the next gen- 
erations of stars in the SC inherit higher metallicities from 
the metal-enriched ISM (Bailin & Harris 2009; Fensch et al. 
2014). Therefore, the more massive a SC is, the higher its 
average stellar metallicity is. Second, readers may note that 
in the Low-metal run, the average metallicity of the DMDG 
is higher than that of the two bound SCs (colored markers in 
the bottom-right panel). A similar discrepancy is observed in 
DF2: [M/H] ~ —1.07 for DF2’s stellar body and ~ —1.63 
for its SCs (Fensch et al. 2019). The higher metallicity in the 
DMDG’s main stellar body could be attributed to the contin- 
uing star formation after the initial starburst. The later gener- 
ations of stars formed in the more metal-enriched ISM raise 
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Figure 6. Top panel: the absolute magnitudes (My) versus the 
half-light radii of the SCs formed in the most massive DMDG in 
the Fiducial run. My is estimated from the SC's mass, assuming 
the mass-to-light ratio M/Ly — 2. The five most massive SCs are 
shown with the colored markers (the same SCs and markers as in 
Figures 2, 4 and 5). We overplot the luminous GCs located in DF2 
(My « —8.5; triangles from van Dokkum et al. 2018c; squares from 
Ma et al. 20202); here the half-light radii and My in the F606W 
band are used. See Section 3.2 for more information. 


the average stellar metallicity of the DMDG, while the SCs’ 
metallicities rarely change once the initial starburst ceases.? 
Finally, we analyze the individual characteristics of the 
SCs, such as their luminosities. Figure 6 shows the scatter 
plot on the plane of the absolute magnitude (My) versus the 
half-light radius (= half-mass radius Rh) for all the SCs in 
the most massive DMDG of the Fiducial run. The abso- 
lute magnitude of a SC is estimated from its mass, assuming 
a typical constant mass-to-light ratio M/Ly = 2 (e.g., van 
Dokkum et al. 2018c). The SC masses plotted in Figure 6 
range from 1.2 x 10? to 9.0 x 106 Mo. Among them, 32 grav- 
itationally bound SCs are with —8.5 « My « —6.5, while 10 
bound SCs are with My « —8.5. For comparison, we in- 
clude the luminous SCs observed in DF2 (My « —8.5; van 
Dokkum et al. 2018c; Ma et al. 2020a). The absolute mag- 
nitudes of the luminous SCs are similar in our Fiducial 
run and DF2. Note that the peak of our cluster luminos- 
ity function (the bottom panel of Figure 6) lies between the 


8 Note that this trend is found only in the Low-metal run, but not in the 


Fiducial run where the initial gas metallicity is already set to 1 Zo. 
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peaks of the Milky Way's GCLF and DF2’s (van Dokkum 
et al. 2018c).? We argue that the extremely dense, shock- 
compressed gas clumps induced by the high-velocity galaxy 
collision could give rise to multiple luminous SCs—a sce- 
nario considered by Trujillo-Gomez et al. (2020, 2021) and 
demonstrated with simulations by Kim et al. (2018) and 
Madau et al. (2020). Meanwhile, the half-light radii Rp of 
the simulated SCs are 1.5-3.0 times smaller than the esti- 
mated radii of the GCs in DF2. As was for Rp of our simu- 
lated DMDG, it is important to note that our simulated SCs 
evolved only for « 800 Myr. If evolved for a sufficiently long 
time, the radii of SCs could change due to, e.g., the dynami- 
cal heating by their host galaxy.!° 


4. CONCLUSION 


The diffuse dwarf galaxies DF2 and DF4, with their al- 
leged absence of dark matter and multiple luminous member 
GCs, have been a topic of intense debate as the astrophysi- 
cal community tries to comprehend their formation mecha- 
nism. Using a suite of high-velocity galaxy collision simu- 
lations with 1.25 pc resolution, we have demonstrated that 
such a collision induces the formation of the DMDGs with 
M, ~ 105 Mz. With numerical resolution significantly im- 
proved from our previous study (Shin et al. 2020), we re- 
solve the dynamical structure of the produced DMDGs and 
find that a co-planar prograde-prograde galaxy collision can 
produce a compact, rotation-supported, oblate DMDG. We 
have also resolved and studied the SCs forming simultane- 
ously with the DMDGs, which later could become progen- 
itors of GCs. A galaxy collision with a high relative ve- 
locity of ~ 300 km s^! creates pockets of severely shock- 
compressed gas that rapidly turn into multiple luminous SCs 
(M, > 106 Mo, My € —8.5 mag) that are in line with the ob- 
served GCs in DF2. The majority of stars in the DMDGs and 
their SCs form less than 150 Myr after the progenitor galax- 
ies’ pericentric approach. Many properties of the observed 


DMDGs’ GCs are reproduced. In our Fiducial run, 10 lu- 
minous, gravitationally bound, virialized SCs have a line-of- 
sight velocity dispersion 11.2 km s™!. In the Low-metal run 
with the initial gas metallicity of 107? Zo, the DMDG shows 
a higher average stellar metallicity than the SCs. 

Our study suggests that a high-velocity galaxy collision 
is a promising candidate for the simultaneous formation of 
DMDGs and their luminous member GCs. In the forthcom- 
ing studies, we aim to examine the galaxy collision-induced 
DMDG formation scenario with more sophisticated subgrid 
physics that fully takes into account the evolution of inter- 
stellar medium due to energetic stellar feedback with even 
higher resolution (see footnote 2), in a cosmological con- 
text for ~ 10 Gyr, comparable to DF2’s age (~ 9 Gyr; van 
Dokkum et al. 2018c; Fensch et al. 2019) designed to follow 
its later evolution. 
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